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Abstract 

The equilibrium of an axisymmetric magnetically confined plasma with anisotropic 
resistivity and toroidal flow is investigated in the framework of magnetohydrodynamics 
(MHD). The stationary states are determined by an elliptic differential equation for the 
poloidal magnetic flux function tp, a Bernoulli equation for the pressure and two relations 
for the resistivity components r)\\ and rj± parallel and perpendicular to the magnetic field. 
The flow can affect the equilibrium properties solely in the presence of toroidicity because 
in the limit of infinite aspect ratio the axial velocity does not appear in the equilibrium 
equations. The equilibrium characteristics of a tokamak with rectangular cross-section 
are studied by means of eigenf unctions in connection with exact solutions for the cases of 
"compressible" flows with constant temperature, T(ip), but varying density on magnetic 
surfaces and incompressible ones with constant density, g(ip), but varying temperature 
thereon. Those eigenf unctions can describe either single or multiple toroidal configura- 
tions. In the former case the equilibrium has the following characteristics: (i) the rj\\- 
and r/_i_-profiles on the poloidal cross-section having a minimum close to the magnetic 
axis, taking large values on the boundary and satisfying the relation rj± > m are roughly 
collisional (ii) the electric field perpendicular to the magnetic surfaces possesses two local 
extrema within the plasma and vanishes on the boundary and (iii) the toroidal current 
density is peaked close to the magnetic axis and vanishes on the boundary. The impact 
of the flow and the aspect ratio on the aforementioned quantities is evaluated for both 
"compressible" and incompressible flows. 
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1. Introduction 



Understanding the equilibrium properties of a magnetically confined plasma is one 
of the key issues in the fusion research. The majority of equilibrium studies up to date 
concern plasmas of zero electrical resistivity and have been based on the Grad-Schluter- 
Shafranov equation which describes the magnetohydrodynamic (MHD) equilibria of ax- 
isymmetric plasmas. Additional physical input in terms of resistivity and plasma flow, 
however, are very important. In addition to its apparent role for Ohmic heating, the im- 
portance of resistivity is connected to the operation of a steady-state fusion reactor which 
will involve time scales much longer than the resistive-MHD one. Also an attractive fea- 
ture of such a reactor would be operation under minimum sources of mass, momentum 
and energy. Sheared flow plays a key role in the formation of both edge transport barriers 
(L-H transition) |2] and internal transport barriers 0-11] i n tokamaks. This flow is 
associated with sheared electric field profiles [7|-PII| with significant amplitudes in the 
barrier region. Another pertinent quantity is the safety factor closely related to 

the toroidal current density. 

Theoretically it was proved long time ago [L4J that axisymmetric resistive MHD equi- 
libria with scalar resistivity uniform on magnetic surfaces are not compatible with the 
Grad-Schluter-Shafranov equation and the poloidal current density should vanish. The 
nonexistence of axisymmetric equilibria with constant resistivity was also suggested in 
Refs. J03 EH! • To examine whether these undesirable properties can be removed by in- 
cluding, in addition to resistivity in Ohm's law, flow and viscosity terms in the momentum 
equation is a formidable task, which should be preceded by a step by step thorough inves- 
tigation of particular cases possibly introducing additional physical input each time. In 
this respect two of the authors studied axisymmetric equilibria with scalar resistivity and 
flow purely toroidal jT7j and parallel to the magnetic field B ^Sj and found that the flows 
considered can not remove the aforementioned inconsistencies. Non-vanishing poloidal 
currents are possible in steady states with parallel flows in the presence of anisotropy, 
i.e. for different resistivity components 7711 and rj± parallel and perpendicular to B 19 ; 
however, in this case neither 7711 nor rjx can be uniform on magnetic surfaces. The sole 
external source in Refs. * s ^ e toroidal current loop voltage. Resistive equilibria 

in a similar spirit of minimum external sources were investigated in Refs. (THJ EI] and |2Uj - 
[2~lj . The particular flow directions considered in Refs. JH-jlH] are no ^ inconsistent with 
Pfisrch-Schluter diffusion in the sense that the equilibrium solutions constructed therein 
neither exclude nor can be included in possible Pfirsch-Schluter-diffusion solutions (hav- 
ing velocity components perpendicular to the magnetic surfaces). Also, it is reminded 
that the pertinent pioneering study [2H] does not include external sources of current and 
neglects the flow term in the momentum equation. 

In this report we extend the studies ^21-^3 to equilibria with toroidal flows and 
anisotropic resistivity. There are two advantageous features of toroidal flows compared 
with parallel ones: (i) they are associated with electric fields E p perpendicular to the 
magnetic surfaces and (ii) as we will show exact solutions with magnetic surfaces of uni- 
form temperature, T = T(ip), are possible. Equilibria with incompressible flows having 
uniform density but varying temperature on magnetic surfaces will also be examined. In 
this respect it is noted that although purely toroidal axisymmetric flows are inherently 
incompressible because of symmetry, T(?/>)-equilibria can be regarded as "compressible" 
in the sense that the density varies on magnetic surfaces. In both cases the study can be 



2 



carried out analytically up to the construction of exact solutions. In particular equilibrium 
eigenstates of a tokamak with rectangular cross-section will be derived in connection with 
exact solutions describing either single toroidal or multitoroidal configurations. Further- 
more for single toroidal eigenstates we will study the characteristics of the conductivity 
components cry and a±, the electric field E p , and the toroidal current density along 
with the impact of the flow on these quantities. This impact will be examined by varying 
a sound-speed Mach number M for "compressible" flows and a parameter A relating 
to the q- and E p - profiles and their variation perpendicular to the magnetic surfaces for 
incompressible ones. In particular the uniformity of cry and a± on magnetic surfaces will 
be examined independently of solutions. Also it will be shown that the impact of the flow 
on the equilibrium is crucially related to the toroidicity. 

An outline of the report is as follows. Reduced equilibrium equations involving if), the 
pressure and the resistivity components are derived in Sec. 2. Tokamak eigenstates are 
constructed in Sec. 3 for both "compressible" and incompressible flows. The equilibrium 
characteristics of a±, an, E p , and along with the impact of the flow and the aspect 
ratio on them is the subject of Sec. 4. The possible role of the flow shear is also briefly 
discussed therein. The conclusions are summarized in Sec. 5. 



2. Equilibrium equations 

In this section we shall derive reduced equilibrium equations for an axisymmetric 
magnetically confined plasma with anisotropic resistivity and toroidal flow. The procedure 
is unified in the sense that relevant energy equations or equations of state are not adopted 
from the beginning; they will specified when necessary later. 

The starting equations in standard notation and convenient units are the following: 

V • (f?v) = 0, (1) 

g (v . V)v = J x B — VP, (2) 

V x E = 0, (3) 

V x B = J, (4) 
V ■ B = 0, (5) 

E + v x B = T) ■ J = 77|| ■ J || + T]±J±, (6) 
An energy equation or equation of state , (7) 

where 

is the resistivity tensor; the indices || and _L indicate directions parallel and perpendicular 
to B; accordingly Jm = (J • b)b, = bx (Jxb) = J — J with b = B/P. The procedure 
to follow is based on identifying some integrals as flux functions, i.e. functions constant 
on magnetic surfaces, and reducing the set of Eqs. (PEJ) by projecting the momentum 
equation (J2J) and Ohm's law © along the toroidal direction, the poloidal one (or parallel to 
B when convenient), and perpendicular to the magnetic surfaces. Important information 
is also drawn from an integral form of © [Eq. (|13|) below]. 

In cylindrical coordinates (R, <f>, z) with z corresponding to the axis of symmetry the 
equilibrium quantities for the case under consideration do not depend on the toroidal 
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angle 0; the toroidal velocity and the divergence-free magnetic field and current density 
can be expressed, with the aid of Ampere's law, in terms of the functions K(R, z), ^(R, z) 
and I(R, z) as: 

£ V = (9) 
B = /V0 + V0xV^, (10) 
J = A*ipV<f) - V(f) x V/, (11) 

where ip labels the magnetic surfaces and A* = R 2 V -(V/i? 2 ). 

By projecting the momentum equation along the toroidal direction one obtains 

V(f) ■ (W x V/) = 0, (12) 

which implies that / = Therefor, unlike the case of parallel flow the current 

surfaces coincide with the magnetic ones irrespective of equation of state. Integration of 
© along a contour c defined by the cut of an arbitrary current surface with the poloidal 
plane yields the equation: 

J~E- d\ + J (v x~B) ■ dl = J(rj ■ J) • dl, (13) 

where dl = V0 x V^/| V0 x is the unit vector along the poloidal direction. Since in 
equilibrium it holds that d~B/dt = 0, the first integral on the left-hand side of ()13|) vanishes 
by Stoke's theorem. Also the second integral vanishes due to the toroidal direction of the 
flow. For the integral on the right-hand side to vanish the integrand must necessarily do 
so because the J po /-lines are closed, nested and V ■ J = 0; therefor it should hold locally 

(77 • J) • dl = (rj ■ J) pol = 0. (14) 

For isotropic resistivity, = 7711, (|T4^) implies that the poloidal current density must 
vanish. In the presence of anisotropy, however, non-vanishing poloidal current densities 
are possible as expected because the toroidal electric field can drive a current in the 
poloidal direction. The rest of the report concerns equilibria with non-vanishing poloidal 
current densities. 

Expressing the electric field on the poloidal cross-section in terms of the electrostatic 
potential, E p = — V$, the component of local Ohm's law © in the poloidal direction on 
account of ((Hj) yields 

V0 ■ (V$ x VV>) = 0. (15) 

which implies that $ = $ (?/>); therefor E p is perpendicular to the magnetic surfaces. The 
total electric field is given by 

E = V c V(f) + E p = V c V(j> - $'VV>, 

where 2ttV c is the constant toroidal loop voltage and the prime denotes derivative with 
respect to ip. Subsequently, the component of © along yields 

($'-^HW| 2 = 0, (16) 
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and therefor the quantity 



K 

^R2 



(17) 



identified as the rotation frequency, is a flux function uj = u)(ip). Eq. (JT4*j) and the 
component of (JHJ) in the toroidal direction respectively yield the following equations: 

Arj 



. -(/A*V-/'|VVf)-^I' = 0, 



(18) 
(19) 



where Arj = r]± — r^y . Any equilibrium solution should be compatible with (fTHj) and ([!§]) 
which, accordingly, can be solved for T]± and 7711 to yield 



Vr 



A*ip + LP 
I'(BR) 2 
IA*ip - J'IV-0 



(20) 
(21) 



With the aid of the integrals 7 = I(ip), $ = $('0) and to = oo(ijj) the components of 
along B and respectively yield 

'uj 2 R 2 - 



VP 



[A'V + I/'JlVVf + R 2 



Q v 2 

2- 



•B = 0, 



= 0. 



22 



(23) 



Owing to the axisymmetry and the toroidal direction of the flow these equations do not 
contain the resistivity and are identical in form with the respective ideal-MHD equations. 

In order to reduce Eqs. (|2*2^1 and PHJ) further an energy equation or equation of state 
is necessary. Owing to the large heat conduction along B, isothermal magnetic surfaces, 
T = T(ip), is an appropriate equation of state for fusion plasmas. In this case employing 
the ideal gas law, P = XgT, integration of (J22j) yields 



P = PsW exp( 



u 2 R 2 
2AT 



where P s {ip) is the pressure in the absence of flow. With the aid of 
to the final "compressible" equation 



Eq. 



A> + II' + R 2 



,+ ' 2 \\t) 



exp 



u 2 R 2 
2XT 



0. 



(24) 
leads 

(25) 



For ideal plasmas (fj—0) Eq. (J25)) was originally obtained in Ref. 
An alternative equation of state is incompressibility: 

V • v = 0. 



(26) 



Consequently, (JTJ) implies that the density is a flux function, g = g(ip), and therefor one 
can find along the same lines the following incompressible equations for P and ip: 

R 2 gu 2 



P = P,W) + 



(27) 
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a*v + //' + r 2 p' s + ~Y^ 2 y = °- ( 2s ) 

Eq. ()28j) is identical with a particular form of the axisymmetric equilibrium equation for 
incompressible flow of arbitrary direction obtained in Ref. j2Z] for ideal plasmas. 

Once Eqs. (|25|) and ()28|) are solved for ip the resistivity components can be determined 
by (j2Ti|) and (j2~Tj) . In general inspection of Eqs. (|2T!|) and (}2*T]) implies, like the case of 
parallel flows [T^], that neither ?7_|_ nor 7711 can be uniform on magnetic surfaces; indeed 
solving (j23|) and (|28|) for A*^ and substituting into (J2T?|) and (|2*T|). one can see that 
and m depend, in addition to ip, explicitly on R (and on |V?/>| 2 as concerns 7711). However 
as we will see in Sec. 4, and 7711 can be collisional-like, viz. they can have a minimum 
close to the magnetic axis, take very large values on the boundary and it holds that 
i]± > m. Reasons for temperature deviations on magnetic surfaces, which can result to 
non uniformity of r}\\ and rjj_ thereon, are discussed in Sec. Ill of Ref. |19j . 

Summarizing this section, the MHD equilibrium states of an axisymmetric plasma 
with anisotropic resistivity and toroidal flow is governed by an elliptic differential equa- 
tion for the poloidal magnetic flux function [Eq. (|2*3|l for "compressible" flow and (|28|) for 
incompressible one], a Bernoulli relation for the pressure and self-consistent expressions 
for the resistivities 7711 and T]x- Both Eqs. (|25|) and (|28|) contain four flux-functions, three 
out of which, i.e. P s , I and lu, are common. The fourth function is T for the "compressible 
equation" and g for the incompressible one. For vanishing flow ()25j) and ()28|) reduce to 
the Grad-Schliiter-Shafranov equation. 



3. Exact solutions 

Linearized forms of Eqs. (|2*3jl and ()28|) in connection with appropriate assignments of 
the free flux functions they contain can be solved analytically. In the present study we 
will employ exact solutions as follows. 

"Compressible" flow 

The ansatz used to linearize Eq. (|23|) is [2E| [2H] 

I 2 = II + Ili) 2 

P s = 2P ^ 2 (29) 

uj 2 7 M 2 

— = = constant 

Here, Iq/R is the toroidal vacuum field, the parameter Ii describes the magnetic properties 
of the plasma; P , 7, and M are a pressure parameter, the ratio of specific heats, and 
the Mach number with respect to the sound-speed at a reference point (z — 0, R — R ) 
with R to be specified later. Note that the toroidal current density profile can vanish on 
the plasma boundary via ([TUjl . 

Eq. (J25J then has a separable solution, lZ(R)Z(z), when the constant of separation is 
equal to RqI\. For configurations symmetric with respect to mid-plane z = this solution 
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is written in the form 



ip(x,y) = d 



Jo 



2t\ e 



7 Mq 2 



7 M 



C2Y0 



2t\ e 



7 M 



cos(i? Az/), 



(30) 



where x = R/Ro and y = z/R ; J and Yq are zeroth-order Bessel functions of first- and 
second-kind respectively; and r 2 = 4P -Ro- 

Incompressible flow 

In this case the ansatz employed to linearize (J28|) is (201 



P = % + I 2 ^ 2 



(31) 



K 2 
qR* 



The third of equations (|3*T|) on account of (|T7|) indicates that A is associated with the 
density and electric field profiles and their variation (shear) perpendicular to the magnetic 
surfaces [{quo 2 ) 1 7^ 0]. The polarity of E p and the aforementioned shear permits A to take 
either positive or negative values. This is a remarkable difference as compared with the 
choice (|2E1) which is shearless [{u 2 /\T)' = 0]. Also, note that, unlike M in ASH, A is 
dimensional. 

A separable solution is now expressed in terms of the first- and second-kind Airy 
functions, Aj and E>i, as [3U] 



ip(x,y) = Ci 



Ai 



PiR t 



cos (Rohy). 



(32) 



In connection with solutions (pffij) and (|H2*|) we are interested in the steady states of a 
tokamak the plasma of which is bounded by a conducting wall of rectangular cross-section, 
as shown in Fig. ^ In addition, we assume that the plasma boundary coincides with the 
outermost magnetic surface. Thus, the magnetic field is tangential to and the pressure P 
must vanish on the boundary [31]; accordingly, the function if) should satisfy the following 
boundary conditions 

V(j/±) = (33) 



and 



ip(x z 



0. 



(34) 



where y± = ±a/Ro and x± = 1 ± bj Rq. The equilibrium becomes then a boundary- value 
problem. Eigenstates can be determined by imposing conditions (|33J) and (|34p directly to 
solutions ({HO} and (|3*2^1 . Specifically, applied to the z-dependent part of the solutions 
yields the eigenvalues 

J ' = a"( £7r "i)' i=1 > 2 >- ( 35 ) 
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Figure 1: The cross-section of the plasma boundary. The aspect ratio a is defined as 
Ro/b, where R is the geometric center of the configuration. 



for the quantity Ji which is related to the poloidal current function The respective 

eigenfunctions are associated with configurations possessing i magnetic axes parallel to 
the axis of symmetry. Condition (J34)) is pertinent to the independent part of the solution. 
Owing to the flow this part contains the parameter M in the "compressible" case and 
A in the incompressible one in addition to the pressure parameter Pq. Thus, condition 
(134)1 can determine flow eigenvalues depending on the parameter Pq which remains free, 
F n (P ) (n = 1,2,3,...) with F standing for either M or A, or vice versa, pressure 
eigenvalues Pq(F) with F being free. The other parameters C\ and Ci contained in ()3U)1 
and (J32)) are adapted to normalize ip with respect to the magnetic axis and to satisfy the 
boundary condition (|34)l respectively. The eigenfunctions in association with F n (Po) [or 
Pq(F)] are connected to configurations having n magnetic axes perpendicular to the axis 
of symmetry. Therefor the total equilibrium eigenfunctions ipin — Z^{z)TZ n {R) describe 
multitoroidal configurations having i x n magnetic axes. 

On the basis of the above solutions one can evaluate the impact of the flow on the 
resistivity components r/± [Eq. (|2U|) ] and r}\\ [Eq. (|2*T|) ]. the electric field perpendicular 
to the magnetic surfaces [E p = — V$ = — $'VV>] and the toroidal current density 
[Eq. ()TT)) ] for both "compressible" and incompressible flows. It is emphasized that this 
impact is crucially related to the toroidicity because in the limit of infinite aspect ratio 
the equilibrium equations do not contain the axial velocity regardless of "compressibility" . 
Indeed, for a cylindrical plasma of arbitrary cross-section the equations respective to (|22)l 
and (|23|l read 

B ■ VP = (36) 
V>+(p+Sy = o. (37) 

For ideal plasmas Eqs. and (pT7)) follow respectively from (16) and (17) of Ref [32| 
for vanishing poloidal velocity {F' = therein). Therefor the flow may have an impact 
on equilibrium only in the presence of toroidicity. Also, note that the pressure becomes a 
flux-function. Because of the importance of toroidicity the impact of the aspect ratio on 
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equilibrium in addition to that of the flow will be evaluated in the next section. 



4. Impact of the flow and aspect ratio on equilibrium 

The quantities to be examined are the conductivity components a± = l/i]± and a\\ = 
I/7711, the electric field E p and the toroidal current density 3$. Both "compressible" and 
incompressible flows will be studied for single toroidal configurations. The study is based 
on the eigenfunction which for "compressible" flow is shown in Fig. |2]and will be made 
by varying the flow parameters Mo and A for "compressible" and incompressible flows 
respectively. For any value of M or A the respective lowest eigenvalue of the pressure 



z 



1 




-1 



Figure 2: Magnetic surface cross-sections for the single toroidal configuration in connec- 
tion with the eigenfunction ipu for M = 0.5 and a = 3. 



parameter P will be calculated numerically. The variation of M and A will correspond to 
the same range of eigenvalues of Pq. Specifically, for aspect ratio a = 3 the flow parameters 
Mo and A will be ranged in the intervals [0.1, 0.7] and [-0.001, -0.01] respectively. For 
a = 2 the respective intervals will be [0.1, 1] and [-0.001, -0.08] unless stated otherwise. 
Variation of the flow parameters in connection with the results to be presented will refer 
to these intervals. The tokamak scaling B p w 0.1-B^, will be used for the calculations. 
Also, we should make the following clarification: it occurs that solutions (jBUj) and (}3*2*j) 
oscillate as the flow parameters are varied, viz. for any given point (R, z) these solutions 
considered as functions ^(M )and ip(A) take successively larger and lower values as M 
and A are varied monotonically. This, would give rise to an oscillatory behavior to all 
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physical quantities which is physically unjustifiable. For the conductivity components 
this can be seen in Fig. El To avoid this difficulty the solutions will be normalized in 




Figure 3: A set of curves demonstrating the oscillation of the profile of a± on the mid- 
plane z — 0, normalized with respect to a constant value a c , for "compressible" flow when 
the value of the function tp on the magnetic axis is flow dependent and the Mach number 
Mq increases: a) Mq = 0.1, b) M = 0.2, c) Mq = 0.5 and d) M = 0.6. The aspect ratio 
is a = 2; x = R/Ro with the vertical axis being placed at the position of the geometric 
center (x = 1). 



such a way that the poloidal magnetic flux on the magnetic axis is unity irrespective of 
flow. This is accomplished by choosing appropriately the parameter C\ for each value 
of Mq or A [Ci(Mq) for "compressible" flow and C\(A) for "incompressible" one]. Also, 
solutions ([30)1 and (|32|1 have a strong parametric dependence on the flow parameters 
[note the exponential dependence of (jSOjl on Mq]. Consequently, this dependence results 
in large quantitative changes in the physical quantities for large M (M ~ 1) or small 
A (A m —0.01) which most probably overestimate the impact of flow. In addition, it is 
noted that except for the conductivity components an increase of Mq has qualitatively the 
same impact on the physical quantities to be examined with that caused by a decrease of 
A. 

The results concerning the characteristics of the quantities o"j_, an, E p , and and the 
impact of the flow and a on them are as follows. 

Conductivity components 

1. For both "compressible" and incompressible flows the profiles of a± and o~\\ on the 
poloidal cross-section are collisional-like, i.e. they have a maximum close to the 
magnetic axis vanish on the boundary and stands in most of the cases that an > 
o~±^M\- Profiles of the resistivity components on the mid-plane z = are shown in 
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Figure 4: The figure shows the increase of the maximum of a± and ay and the displacement 
of its position outwards for "compressible" flow due to the increase of the Mach-number 
M . 




Fig. H]for "compressible" flows and in Fig. EJfor incompressible ones. For vanishing 
flow and a = 3 one finds Ac/cry = (cry — <t±)/<t\\ = 0.2122; for the reversed-field- 
pinch scaling, ~ B p , the value of (cry — cr±)/(T\\ becomes double, i.e. 0.4244. Also, 
for a = 3 increase of M from 0.1 to 0.7 results in a percentage decrease of Act by 
4% while decrease of A from -0.001 to -0.01 leads to an increase of Act by 3.4%. 

2. The maximum of the profiles takes larger values as Mq increases and lower values as 
A decreases. For a = 3 increase of Mq (from 0.1 to 0.7) leads to percentage increases 
of cr_i_ and cry by 9% and 7% respectively. For incompressible flows the respective 
decreases are 73% and 67% (for a = 3 and decrease of A from -0.001 to -0.01). 

3. The position of the maximum is shifted outwards from the axis of symmetry as Mq 
increases and is nearly insensitive to the variation of A [32] (see Figs. 0] and 03). For 
example for a = 3 the position of the maximum of <jj_ changes from 1.119 to 1.151 
(when Mq varies from 0.1 to 0.7). 
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4. The lower a 

• the larger the decrease of Act = mi — a± as M increases; 

• the smaller the increase of Act as A decreases; 

• the larger the increase of the a±- and ay-maximum as M increases; 

• the smaller the decrease of the a±- and ay-maximum as A decreases; 

• the larger the shift of the position of a±- and ay-maximum as Mq increases. 

Those conclusions become evident if the following results for a = 2 are compared 
with those presented above for a = 3: 

• the decrease of Aa for compressible flow is 13%; 

• the increase of Aa for incompressible flow is 1%. The respective variation of 
A in this case is [-0.001,-0.01] (see [M]); 

• the increase of the conductivity maximum values for "compressible" flows are 
44% (a ± ) and 27% (a,,); 

• the decrease of the conductivity maximum values for incompressible flows are 
43.7% (a±) and 63% (ay); 

• the position of the maximum of a± for "compressible" flow is displaced from 
1.228 to 1.33. 

The "compressibility" dependent impact of the flow on a±, an and Aa may relate to the 
flow shear of the incompressible solution (|32|) while the "compressible" solution (|3U|) is 
shearless. 

Electric field 

Prescribing the rotation frequency by 

u = co ip n , (38) 

where n is a shaping parameter, the electric field, E p = — $'V-?/>, with the help of (J 17)) 
becomes 

E p = -to ip n Vip. (39) 

The profile of E p on the poloidal cross-section possesses two local extrema within the 
plasma volume one on the left-hand-side and the other on the right-hand-side of the 
magnetic axis with opposite sign, which means that the electric filed has different polarity 
[as can also be deduced from Eq. ([39p] and is similar to ones observed in equilibria with 
internal transport barriers. The profile \E r \, on the mid-plane y = is shown in Fig. 
Increase of M Q or decrease of A leads to an increase of both local maximum values of \E r \ 
with the one outwards more than the one inwards the magnetic axis and to a displacement 
of their positions outwards. For "compressible" flows this is shown in (Fig. IHJ). For a = 3 
the maximum values of \E r \ increase by 4% for "compressible" flow (increase of Mq from 
0.1 to 0.7) and 5.6% for incompressible one (decrease of A from -0.001 to -0.01). The 
lower a the larger the increase of the maxima and the displacement of their positions 
outwards. As an example, for a = 2 the increase of the maximum \E r \ becomes as large 
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X 



0.6 0.8 1.2 1.4 



Figure 6: Electric field profiles on the mid-plane z = for "compressible" flow Comparison 
of the two curves indicates the increase of the local maxima of E p and the outward shift of 
the their positions as the Mach number Mo increases from 0.8 to 1. The point in between 
the two maxima at which E p = corresponds to the magnetic axis. The aspect ratio is 
a = 2 and the value of the shaping parameter n [Eq. (|58|)] is 3. 



as 15% for "compressible" flow with strong increase occurring above Mq = 0.8 and 9% 
for incompressible flow. Also, when the parameter n takes larger values the maxima of 
\E r \ take lower values but the profile of E r becomes steeper and more localized (Fig. |7j); 
therefor the shear Se t = dE r /dx is also increased. 




x 



0.6 0.8 1.2 



Figure 7: Electric field profiles on the mid-plane z = for a = 3, Mq = 0.4 and two 
different values of the shaping parameter n: n = 1 and n = 3. 

The toroidal current density 
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The profile of the toroidal current density, 

J = ^A>, (40) 

is peaked with its maximum in the vicinity of the magnetic axis and vanishes on the 
plasma boundary (Fig. [SJ). Like E p , as Mq increases or A decreases the maximum of 

(J*)c 

1.2 b _ 




0.6 0.8 1.2 1.4 

Figure 8: Toroidal current density profiles on the mid-plane z = for "compressible" flows 
Comparison of the two curves indicates the increase of the maximum and the outward 
displacement of its position as the Mach number Mq increases from 0.1 to 0.8. The aspect 
ratio is a = 2. 

takes larger values and its position is displaced outwards. Also, the impact of the varia- 
tion of the flow parameters is stronger as the aspect ratio becomes smaller. In particular 
for a = 3 and a = 2 the maximum value of the current density increases by 8% and 
35.5% respectively for compressible flows (increase of Mq). The respective increases for 
incompressible flows (decrease of A) are 5.6% and 18%. 

5. Conclusions 

In this report we have investigated the MHD equilibrium of an axisymmetric mag- 
netically confined plasma with anisotropic resistivity and toroidal flow subjected to the 
single external source of toroidal current loop voltage by including the flow term in the 
momentum equation. Equilibria of this kind are inherently free of Pfirsch-Schliiter diffu- 
sion. Also, unlike the case of parallel flows, there is an electric field E p perpendicular to 
the magnetic surfaces. The study includes "compressible" flows with varying density but 
uniform temperature on magnetic surfaces and incompressible ones with uniform density 
but varying temperature thereon. It turns out that the equilibrium states are determined 
by an elliptic differential equation for the poloidal magnetic flux function ip [Eq. (|25|) for 
compressible flows and (J28j) for incompressible one], a Bernoulli equation for the pressure 
[respective Eqs. (J2U) and (|27j) ] and two relations for the resistivities r)± [Eq. (|2*UJl ] and 
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r)\\ [Eq. (|2ip] parallel and perpendicular to the magnetic field. Owing to axisymmetry 
and the toroidal direction of the flow, the equilibrium equations and pressure relations 
are identical in form with the respective ideal-MHD ones. The impact of the flow on 
equilibrium can be "activated" solely in the presence of toroidicity because the cylindrical 
equilibrium equations do not contain the axial velocity. 

The equilibrium of a tokamak plasma bounded by a rectangular cross-section has 
been studied by means of equilibrium eigenfunctions in connection with exact solutions 
for "compressible" and incompressible flows. These eigenfunctions can describe either 
single toroidal or multiple toroidal configurations. In the former case we have studied 
the characteristics of the conductivities a± and cry, the electric field E p , and the toroidal 
current density as well as how they are affected by the flow. The impact of the flow has 
been examined by varying pertinent flow parameters, i.e. a sound-speed Mach number 
M Q for the "compressible" flow and a parameter A relating to the density and electric 
field and their variations perpendicular to the magnetic surfaces for the incompressible 
one. For (i) single toroidal configurations (ii) eigenfunctions normalized so that ip is 
unity on the magnetic axis regardless of flow (a normalization made to avoid physically 
unjustifiable oscillation of the solutions on the flow parameters) and (iii) variation of 
the flow parameters corresponding to the same variation of the lowest eigenvalue for the 
pressure parameter Pq we came to the following conclusions: 

1. For a toroidal frequency peaked on the magnetic axis and vanishing on the 
boundary [Eq. (JHHj) ] the profile of |E P | on the poloidal cross-section has two maxima 
located the one on the left-hand-side and the other on the right-hand-side of the 
magnetic axis and it vanishes on the boundary. When the maximum of u becomes 
larger and its profile more localized, the profile of |E p | becomes as well more localized 
though its maxima become smaller. 

2. The profile of is peaked with its maximum close to the magnetic axis and vanishes 
on the boundary. 

3. Although the conductivity components can not be uniform on magnetic surfaces 
(this follows in general by inspection of the equilibrium equations) their profiles are 
roughly collisional, viz. they have a maximum in the vicinity of the magnetic axis, 
vanish on the boundary and it holds, in most of the cases that cry > a±. 

4. As M Q increases or A decreases the local maxima of E p and take larger values 
and their positions are shifted outwards from the axis of symmetry. 

5. The impact of the variation of the flow parameters on a± and an rely on "compress- 
ibility": like E p and J^, the maxima of a± and ay become larger as M increases 
and their positions are displaced outwards but these maxima become smaller and 
their positions are nearly not affected as A decreases. Also, the larger the M the 
smaller a\\ — a± but the smaller A the larger a\\ — a±. 

6. For a given value of Mo, the lower the aspect ratio a the smaller the maxima of a\\, 
<t_i_ and Jfj, but the larger the maximum of E p . 

7. For a given value of A, the lower the aspect ratio a the larger the maxima of an , a± 
and E p but the smaller the maximum of J^. 
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8. For increase of Mq or decrease of A (corresponding to the same variation of the 
lowest eigenvalue of Pq) , the lower a the higher the variation of the maxima of E p 
and Jtf, and the displacements of their positions outwards. The impact of a on the 
conductivities are "compressibility" dependent: the smaller the a (i) the larger the 
variation of the maximum of a±, au and cry — <j± (when M Q increases) but (ii) the 
smaller the maximum of a±, a» and au — a± (when A decreases). 

Qualitatively, except for the conductivity components the impact of the flow on E p 
and are independent of "compressibility". The dependence of the results for a± and 
a\\ on "compressibility" may be due to the fact that the incompressible solution (}3*2"j) has 
finite flow shear while the flow of the compressible one (JHUj) is shearless. Quantitatively, 
for a = 2 and increase of Mo from 0.1 to 0.5 or decrease of A from -0.001 to -0.006 
result in percentage variations of all quantities (cry, a±, E p , and J^) less that 10%. Larger 
variations of these quantities associated with higher values of Mq or lower values of A, we 
have found on the basis of solutions ()30|) and ()32|). most probably overestimate the actual 
impact of the flow. 
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